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Soliton interactions of the Kadomtsev-Petviashvili 
equation and generation of large-amplitude water waves 

By Gino Biondini^, Ken-ichi Maruno^ , Masayuki Oikawa^ and Hidekazu Tsujfi 



C/^ ' We study the maximum wave amplitude produced by line-soliton interactions of the 

Kadomtsev-Petviashvili II (KPII) equation, and we discuss a mechanism of genera- 
tion of large amplitude shallow water waves by multi-soliton interactions of KPII. We 

m 

' also describe a method to predict the possible maximum wave amplitude from asymp- 

^ I totic data. Finally, we report on numerical simulations of multi-soliton complexes of 

vn , 

^ I the KPII equation which verify the robustness of all types of soliton interactions and 

o , 

G> , web-hke structure. 

o 

> 

^ ' 1. Introduction 

The existence of waves of large height on the sea surface is a dangerous phenomenon 
inn [T6l I20I l26l [32I . Extreme waves occur much more frequently than it might be 
expected from surface wave statistics |[T6l . These extreme waves, which are particu- 
larly steep and may arise both in deep water and in shallow water, have a significant 
impact on the safety of people and infrastructure, and are responsible for the ero- 
sion of coastlines and sea bottoms and changes to the biological environment. Thus, 
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understanding the physics of these extreme waves is an important task which may 
even contribute to save hves. Ahhough several physical mechanisms of generation 
of extreme waves in deep water have been studied, less is known for the situation of 
shallow water ll26l . 

Waves in shallow water have been studied since the nineteenth century, of course. 
It is well-known that, in the case of weak nonlinearity, for weakly two-dimensional 
cases (that is, for the cases that the scale of variation in the direction normal to 
the propagation direction is much longer than that in the propagation direction) the 
fundamental equations for the dynamics in shallow water may be reduced to the 
Kadomtsev-Petviashvili II (KPII) equation H] [I4l [3T1- Since the KP equation is 
integrable, the theory of integrable system can be used to analyze wave dynamics 
in detail. In particular, Someere et al. recently studied the amplitude of ordinary 
2-soliton solutions of the KPII equation ll27l [33l l34l . They pointed out that the inter- 
action of two solitary waves may be one mechanism of generation of extreme waves. 
In the case of ordinary 2-soliton solutions, the maximum interaction height can be 
four times that of the incoming solitary waves in the limit of the resonant Y-shaped 
solution originally found by Miles Il23ll24ll25l . The interaction pattern of the Miles 
solution and ordinary 2-soliton solutions, however, is stationary. Therefore, it cannot 
describe how a large amplitude wave can be generated. Recently, we studied a class of 
exact line soliton solutions of the KPII equation and found non-stationary interaction 
patterns and web-like structures [5 ]. Such solutions describe resonant line-soliton in- 
teractions which are generalization of the Y-shaped resonant soliton solution. More 
general solutions were later characterized in EllBlffllTlfTSll. 

Our motivation in this paper is severalfold: (i) we determine the maximum ampli- 
tude in interactions of hne sohtons of the KPII equation, (ii) we propose a mechanism 
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to explain the generation of waves of large height, (iii) we describe an algorithm to de- 
termine the maximum amplitude from experimental data, (iv) we perform numerical 
simulations of multi-soliton interactions of the KPII equation, whose results confirm 
the robustness and stability of soliton interactions and web-like structure. 

2. The KPII equation and its soliton solutions 

Here we briefly review some essential results the KPII equation 

{-4Ut + 6UUjc + Uxxx)x + 'iUyy =0 (2. 1) 

that will be referred to through the rest of this work. Hereafter, u{x,y,t) represents 
the dimensionless wave height to leading order, and subscripts x,y,t denote partial 
differentiation. It is well known that solutions of ( 12.11 ) can be expressed via a tau 
function x{x,y,t) as 

u{x,y,t) = 2—\ogT{x,y,t), (2.2) 

where T{x,y,t) satisfies Hirota's bilinear equation lfT2l [29l . Solutions of Hirota's 
equation can be written in terms of the Wronskian determinant ^ [lOl 

t(x,3;,0 =Wr(/i,--. ,/;v) =det(/i"'-'*)i<„y<;v, (2.3) 

with fn^ = djfn/dxj, and where /i , . . . are a set of linearly independent solutions 
of the linear system 

dy dx^ ' dt dx^ 
For example, ordinary A^-soliton solutions are obtained by taking /„ = e^" ' -\-e^^" for 
n = \,...,N, where 

d,„{x,y,t) = -kmX + kl^y - k^J + dmfl (2.5) 
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for m = 1 , . . . , IN, where the 4A'^ parameters ki<-- - < k2N and ^j-q, . . . , QiN-fi are real 
constants. For N =\, one obtains the single-sohton solution of KPII: 

utj{x,y,t) = ^a^^.sech^ [1(0; - Bj)] , (2.6) 

where / = 1 and j = 2. Equation ( 12.61 ) is a traveling-wave solution, [with wavenumber 
k = {kj — ki,kj — kj) and frequency (O = kf — kj], exponentially localized along the 
line 6, = dj of the xy-plane, and is therefore referred to as a line soliton. We refer to 

respectively as the soliton amplitude and direction. Note that dij = tana,^-, where 
is the angle made by the soliton with the positive y-axis (counted clockwise), since 
di — dj = {ki — kj) [—X + dijy — {kf + kikj + kj)t] + dj-o — Oj-q. The actual maximum 
of the solution, or wave height, is given by Uij = jofj. 

The case with = 1 and / = e^' H \-e^", with 0i , . . . , 0m still given by (1231 ). 

was studied in Il22l . In particular M = 3 yields the Y-shaped solution often called 
Miles resonance, in which three line solitons interact at a vertex, as shown in Fig. [T] 
and whose wavenumbers and phase parameters satisfy the resonance conditions ki + 
k2 = k3 and (Oi + (02 = (O3. (Note that, while the Miles resonance is also a travel- 
ing wave solutions, solutions with M > 4 are not ll22l .) More in general, choosing 
fn = f^"^^^ for « = 1 , . . . jA/^ yields T{x,y,t) in the form of a Hankel determinant. In Q 

we studied such solutions with f = e^^ -\ h e^" , and we showed that they produce 

non-stationary and fully resonant (A'^_ jA'^+j -soliton solutions of KPII, that is, solutions 
with = M — N solitons asymptotically as y —00 and = N solitons asymp- 
totically as y ^ 00, and for intermediate values of y these solitons interact resonantly, 
i.e., via fundamental Miles resonances. 
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It should be clear that even more general solutions exist, however. The most 
general linear combination of exponentials can be written as 

M 

fn=Y.(^n,,ne^"\ (2.8) 

m=l 

with 6,„ given by ( 12.51 ) as before. Then ( 12.31 ) is expanded as a sum of exponentials: 

T:{x,y,t) = det{C@K) = Vnn,...MNCmu...,mN exp 0mi,-,m,v ' (2-9) 

l<mi<--<miv <W 

where C = (c„,„,) is the N x M coefficient matrix, = diag(e^' , • • • ,e^"), and the 
M xN matrix K is given by K = {k"^^)- Hereafter, dm,.- -,mN denotes the phase com- 
bination ■■,m„(-x,3',0 = Qmi H 1" ' wMlc V'„!i,...,mA. is the Vandcrmondc de- 
terminant Vmi^...,„,„ =ni <j<j'<N{kmji ~ Kuj) > ^nd C'mi,...,mjv ^hc A'^ X A'^-minor whose 
n-th column is respectively given by the m„-th column of the coefficient matrix for 
n = I, . . . ,N. The only time dependence in the tau function comes from the expo- 
nential phases Qmi....,mN- Also, for all G € GL(A'^,M), the coefficient matrices C and 
C' = GC produce the same solution of KPII. Thus without loss of generality one can 
consider C to be in row-reduced echelon form (RREF), which we will do throughout 
this work. One can also multiply each column of C by an arbitrary positive constant 
which can be absorbed in the definition of Oj-q, . . . , 6m;0- 
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Real nonsingular (positive) solutions of KPII are obtained if ki < ■■ ■ < and 
all minors of C are nonnegative. In f3\ we showed that, under these assumptions and 
some fairly general irreducibility conditions on the coefficient matrix, ( |2.9b produces 
(A^_,A'^+)-soUton solutions of KPII with N-=M -N and A'^+ = N, as in the simpler 
case of fully resonant solutions. Asymptotic line solitons are given by (12.61 ) with the 
indices / and j labeling the phases 0, and Qj being swapped in the transition between 
two dominant phase combinations along the line 0, = Qj. Asymptotic solitons can 
thus be uniquely characterized by an index pair [/, j] with I < i < j < M. We call 
outgoing line solitons those asymptotic as j ^ oo and incoming line solitons those 
asymptotic as j — > — oo. The index pairs are uniquely identified by appropriate rank 
conditions on the minors of the coefficient matrix [3|. We also showed in [31 that 
this decomposition is time-independent, and that the outgoing solitons are identified 
by pairs [i^ ,]„], « = !,... ,N, where i^ .. ,if^ label the A'^ pivot columns of C; sim- 
ilarly, the incoming solitons are identified by pairs [i^ Jn]' " = I, . . . ,M — N, where 
' • ■ • ' Jm-n l^bel the M — N non-pivots columns of C. 

In general, these solutions exhibit a mixture of resonant and non-resonant inter- 
action patterns. In the special case M = 2N, leading to = = N, we call elastic 
A^-soliton solutions those for which the amplitudes and directions of the incoming 
solitons coincide with those of the outgoing solitons, while all other A/^-soliton solu- 
tions are referred to as inelastic. Elastic A/^-soliton solutions were characterized in 
ElllH. The fully resonant line soliton solutions studied in tSJ are a special case of 
the tau function (12.91 ) in which allN xN minors of the coefficient matrix are nonzero. 
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Figure 1: Left: Miles resonance, with {ki,k2,ko,) = (—1,0, 1.5) and f = 0. Center: ordinary 2- 
soliton solution, with (fei , . . . , ^4) = (— 1 , — 0.001 , 0, 1 ) and f = 0. Right: asymmetric 2-soliton 
solution, with {ki,... ,k4) — (—2, —1.5, 1,2) and f = 0. In all cases 0,„;o = form = 1, . . . ,M. 




20 20 



Figure 2: Time evolution of a resonant 2-soliton solution, with (fei, . . . ,^14) = (—2,-1, 1,2) 
and 0„,;o = for m = 1 , ... ,4. Left: t ~ —3; center: f = 0; right: t — 2. 

3. Interaction amplitudes of line-soliton solutions 

Here we first study the amplitude generated by 2-soIitoii interactions. Tliese are ob- 
tained for A'^ = 2 and M = 4, yielding the tau function 

x{x,y,t) = yi,2Ci,2^'^''^ + VlsCls^-^'-^ +yi,4Ci,4^'^'-^ +V23Ci,3e^-' 

+ V2,ACi,Ae^-' + V^,AC^^Ae^'-\ (3.1) 

where now Vij = kj — ki, and with < • • • < as before. The 2-soIiton solutions of 
KPII were recently studied in several papers ll6l[T8l[3ll4ir7l[8]|. 
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3.1. Elastic 2-soliton interactions 



It was shown in 1181 that elastic 2-soliton solutions are classified into three classes, 
shown in Figs. [T]and|2l ordinary (0-type), asymmetric (P-type) and resonant (T-type). 
The coefficient matrices corresponding to these classes have the following RREFs: 

^ \ -\\ 
y 1 1 




Co — , Cp 



Ct 




with C2,3 > C2,4 > 0. These three types of solutions cover disjoint sectors of the 2- 
soliton parameter space of amplitudes and directions 13. Moreover, their interaction 
properties are also different. This difference is obvious in the case of resonant solu- 
tions, but also applies to ordinary and asymmetric solutions, since asymmetric solu- 
tions only exist for unequal amplitude, and the interaction phase shift has the opposite 
sign for ordinary and asymmetric solutions lUl. We next show that these solutions also 
differ in terms of the interaction amplitudes. 

Ordinary 2-solitoii interactions. Consider the coefficient matrix Co in (13.21 ). It is 
Ci^2 = C3.4 = in (13.11 ). while all other minors are unity. In this case, the dominant 
phase combination as x ^ —00 is (2,4) [meaning that e^^-* is the dominant exponential 
in (13. IM . while that as x ^ 00 is (1,3). The asymptotic solitons are [1,2] and [3,4]. 
The pattern of this solution is stationary. The maximum height was discussed in 
Il27ll33l[34l l8l. It is expected that the value of u{x,y,t) at the interaction center would 
give the global maximum. In fact, it was verified in f33l. According to lEl, the value 
of M(x,j,f) at the interaction center, Mq , is given by 

"o''^ = ^(«i,2 + «3,4)+«i,2fl'3,4tanh^^^ , (3.3) 



where 

. , {h-h){k4r-k2) , {d3,A-di;if -{a^,4-ax 2f- ,^ 
Ao = \og- — — — =log— ' • (3.4) 

(^3 - ACl) (k4 - h ) ('33,4 - «1,2 Y " («^3,4 + aij) 

Since (^3 - ^i)(^4 - ^2) - (^3 - ^2)(^4 - ^1) = (^2 - ^i)(^4 - ^3) > 0, > 0. So, 

^(fll 2+^3,4) < "o"^ < ^(«l,2+«3,4)^ (3.5) 

The left hand side is the infinitesimal limit k2 — k\ — > 0+ or /C4 — ^3 ^ 0+ and is 
the sum of the heights of the asymptotic solitons Uy2 and 6^3 4. Accordingly, the 

(ic) 

maximum height Mmax = Mq i^ always greater than the sum of the heights of the 
asymptotic solitons. It is also lesser than Mq"^' := \{aY2 +0^3.4)^ as long as k\ < 
^2 < ^3 < ^4- "max takcs the critical value Uq^'^ in the limit kj,—k2 — > 0+. From 
ai,2fl'3,4 < (flj 2 + 4) /2, for ai,2 = fl'3_4 =: a, Mq'' takes its maximum 2a^ = 4(a^/2) 
which is four times the height of the incoming sohtons. However, it should be noted 
that in the limit kT,—k2^ 0+ the solution degenerates to a Y-shape resonant solution. 

It is also important to note that since 2('^i,2 + '^3,4)^ < 4(max{?7i 2, ?/3^4}) as 
pointed out in [8], the maximum height of u{x,y,t) is less than four times the height 
of the taller(tallest) asymptotic soliton. 

The right hand side of ( 13.51 ) may be also obtained in the following way. Near the 
critical angle (^2 ~ ^3). the interaction arm appears in the intermediate region of 2 
soliton interaction. The interaction arm corresponds to the double-phase transition 
(1,3) ^ (2,4) among the dominant phases in the x3'-plane, which is located along 
the line 0i,3 = 02,4- This is due to the fact that the [l,2]-soliton in 3^ —s- 00 is shifted 
to the right relative to that my ^ —00, while the [3,4]-soliton in 3^ ^ 00 is shifted to 
the left relative to that in y —00. In a neighborhood of the line corresponding to 
the double-phase transition (/i,/2) ^ (71,72). the tau function and the corresponding 
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solution of KP are approximately, up to exponentially smaller terms, 



(3.6) 



and 



(3.7) 



where A,- 



^og{Vj,jJVi,j,) and 



(3.8) 



The double-soliton (13.71 ) is not in itself an exact solution of the KPII equation (and 
hence it is not a true line soliton), because its wavevector and frequency do not satisfy 
the soliton dispersion relation |[T3l . This relation is satisfied in the limit kT,—k2^ 
0+. In this limit there are only three phases in the tau function, and the solution 
degenerates to a Y-j unction. In this limit, the interaction arm of ordinary 2-soliton 
solutions is given by u{x,y,t) = uu24{x,y,t), au24 = a\.2 + fl'3,4- This result agrees 
with the right hand side of ( I3.5I I. This method may be, therefore, effective for deriving 
the supremum for the height of the ordinary 2-soliton solution. 

Asymmetric 2-soliton interactions. Consider now Cp in (13.21 ). It is Ci,4 = C23 = 
while aU other minors are unity. The asymptotic solitons are [1,4] and [2,3], and the 
dominant phase combinations as a; ^ —00 and ;c ^ 00 are (3 , 4) and (1,2), respectively. 
According to lEl, the value of u{x,y,t) at the interaction center, Up'\ is given by 




(3.9) 



where 



4p = log 



(^3 -fcl)(^4 -^2) 
(^2-^0(^4-^3) 



log 



(^1,4 + ^2,3)^ - (^1,4-^2.3)^ 

(ai,4 - ^2,3)^ - (^1,4 - d2,?,Y 




(3.10) 
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Accordingly, we have the estimate 

^(^1,4 -^2,3)^ < 4"''' ^ ^'^1,4- (3-11) 

The left hand side is the resonance limit kj — ki^ 0+ or — — > 0+. The right hand 
side is the infinitesimal limit —k2 ^0^. Note that in this limit a2.3 — ^ 0. It is also 
noted that more accurate estimate for right hand side can be obtained (see HI). Unlike 

(ic) 

the ordinary 2-soliton solution, Mp does not give the global maximum of u{x,y,t). It 
is rather given by the height of the highest asymptotic soli ton ?/i,4 = ^af 4. 

The left hand side of (13.111) can be also obtained in the same way as in the ordinary 
2-soliton solution. For clarity, let us assume ^1 +^4 > /:2 + ^3- In this case, [1,4]- 
soliton is located on the right side of [2,3]-soliton in 3^ — > 00 as in Fig{T] The [1,4]- 
soliton in 3^ ^ 00 is shifted to the right relative to that in 3^ ^ —00, while the [2,3]- 
soliton in J —> 00 is shifted to the left relative to that in y —00. Accordingly, the 
interaction arm appears as the boundary between the two dominant phase regions 
(1,3) iny —00 and (2,4) in y — > 00. Thus, the amplitude of the interaction arm in 
the resonance limit k2 — k\—f 0+ is given by £213.24 as for ordinary solutions. Since 
«13,24 = ^2 + ^4 — (^1 +^3) =k4 — k\ — {ks — ^2) = <^iA — <^2,3, this rcsult agrees with 
the left hand side of (13.111) . For ^1 + ^4 < ^2 + ks, we can also the same final result. 

For ki + k4 = k2 + kj,, [1,4] and [2,3] solitons are parallel and one-dimensional 
over-taking interaction takes place. Even for this case, ( 13.91 ) is still valid if we take 
the interaction center as that in the (x — {ki +k4)y,t)-p\am. According to (13.91 ). the 
height at the interaction center is J7i 4 — 1/2,3- 

Resonant 2-soliton interactions. For the coefficient matrix Cj in (13.21 ) all minors 
are nonzero. The dominant phase combinations as — > —00 and x — > 00 are respec- 
tively (3,4) and (1,2), as with asymmetric solutions, but the asymptotic solitons 
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here are [1,3] and [2,4]. Here the interaction is mediated by four interaction seg- 
ments and the solution is non-stationary. (Its time evolution is shown in Fig. |2l) 
Hence, the situation appears at first to be more complicated than in the other two 
cases. Nonetheless, the calculations are actually simpler, because all the intermedi- 
ate arms are true line solitons of the KPII equation, and the resonance condition is 
satisfied at all vertices. Moreover, the interaction pattern obeys the reflection sym- 
metry {x,y,t) I— > {—X, —y, —t). Indeed, it would be trivial (however, see the last part 
of this subsection) to see that, at all times, the tallest intermediate soli ton is [1,4], 
which is obviously taller than either of the asymptotic solitons. Note, however, that 
ai,3 +^2A — ^^1,4 + '^2,3 > Thus the maximum interaction amplitude is always 
less than the sum of the amplitudes of the asymptotic solitons. In particular, in the 
case of equal-amplitude solitons, ai 3 = (32,4 ='■ fl', we have U^a.^ = 4 < 2a^. Thus 
the maximum height is less than four times that of the asymptotic solitons. 

Figure 2 shows that the interaction at f = neither generate the intermediate soli- 
tons nor large amplitude. At least, we can prove m(0,0,0) < \a\^ for u(x,y,t) with 
zero phase constants dmfl = 0,m = 1, • • • ,4. Our method is as follows: We express 
M (0,0,0), which is a function of kj, in terms of di j and a,j-. It turns out that the terms 
containing di^ are in the form (<i,j- — d,„ nYaijam,n- Accordingly, all the terms can be 
expressed in terms of aij. If we make decomposition like ^1,4 = ^1,2 + ^2,3 +^3,4> we 
find that 4 — m(0,0,0) is a sum of positive terms. 
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3.2. Inelastic 2-soliton interactions 

Inelastic 2-soliton solutions fall into four categories Q, identified by the following 
coefficient matrices in RREF: 



Ci 



Qii 




Ct, 



Civ 




(3.12) 



with r > 0. We consider each of these in turn. For simplicity in what follows we set 
r = 1. Note that in all of these cases exactly one of the minors of C is zero while all 
other minors are unity, which produces a tau function with five phase combinations. 

Type I. In this case Ci,2 = 0. The incoming solitons are [1,2] and [2,4], the outgoing 
solitons are [1,3] and [3,4]. The dominant phase combination as — > — oo is (3,4), 
while that as ;c ^ oo is (1,3). The interaction pattern is a combination of two Y-shape 
resonances, as shown in Figs. [3] and HI In particular. Fig. |4] shows contour plots of 
the solution as well as the index pairs corresponding to each of the intermediate and 
asymptotic solitons. The interaction vertices, however, are not invariant: For f < the 
three phases appearing at each interaction vertex are respectively (2,3,4) and (1,2,3), 
corresponding respectively to solitons [2,3], [3,4], [2,4] and [1,2], [2,3], [1,3]. But at 
f = two resonant stems collide, and the arrangement of solitons changes thereafter, 
with the resonant vertices being characterized by the phases (1,2,4) and (1,3,4) for 
f > 0, corresponding respectively to solitons [1,2], [1,4], [2,4] and [1,3], [1,4], [3,4]. 
(An additional X-shape vertex is produced for f < by the asymptotic solitons [1,2] 
and [3,4]. This interaction is locally the ordinary 2-soliton interaction and the local 
maximum is less than \{ai;i +a^^^Y < ^a\^. We can also prove m(0,0,0) < ^af 4 
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similarly to the case of resonant 2-soliton interactions.) 

The rearrangement in the soliton configuration that happens at f = corresponds 
to the generation of a large-amplitude wave for t > 0. The interaction arm is the 
intermediate soliton [2,3] for ? < and [1,4] for t > 0. The first of these is always 
shorter than the asymptotic solitons [1,3] and [2,4]. The second one, however, is 
taller than any of the others. Moreover, the height of the interaction arm [1,4] is 
f^max = 2ai^4 = ^(fl'ij +ajA)^ > ^(^ij + ' j = 2,3. Thus, the height of the 
soliton [1,4] is always greater than the sum of the heights of the incoming and the 
outgoing solitons. As before, the maximum value of the interaction height relative 
to the height of the asymptotic solitons occurs in the case of equal amplitudes: when 
«i.2 = '32,4 =• a:, it is ?7max = 2fl'^, yielding a ratio of four to one. However, it should 
be noted that since ^ > a, Umax is less than four times the height of [l,3]-soliton. 
That is, the maximum height U^ax is less than four times the height of the highest 
asymptotic soliton. If we impose still further the condition ai 3 = £23,4 = a, the solution 
degenerates to a Y-shape resonant solution. 

Type II. In this case C3,4 = 0. The incoming solitons are [1,3] and [3,4], the outgo- 
ing solitons are [1,2] and [2,4]. The whole solution is a time reversal of the inelastic 
2-soliton solution of type I: that is, uii{x,y,t) = ui(—x, —y, —t). The dominant phase 
combination as ;c — > —00 and x ^ 00 are now respectively (2,4) and (1,2). The tall 
interaction arm is still the soliton [1,4], but now it appears for f < 0. 

Type III. In this case €2,3 = 0. The incoming solitons are [1,4] and [2,3], the outgo- 
ing solitons are [1,3] and [2,4]. The dominant phase combination asx ^ —00 is (3,4), 
while that as x —s- 00 is (1,2). The interaction dynamics is similar to that of solutions 
of type I, as shown in Figs.[5]and[6l The interaction pattern is a combination of two Y- 



14 




Figure 3: Time evolution of the inelastic 2-soliton solution obtained from the coefficient 
matrix d in ( 13.121 ) with (^i, . . . ,k4) = (-2,0,2,4). Left; t = -2; center: t ^ 0; right: t ^ 2. 
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Figure 4: Contour plots of the inelastic 2-soliton solution shown in Fig. |3] The dominant 
phase combinations (oti,OT2) and the index pairs that uniquely identify the line solitons 
are labeled. Left: t — ~\; center: f 0; right: f = L 

shape resonances. The pattern differs depending on the magnitude relation between 
k\ + kn and ^2 + ^^3- Figures |5] and [6] show the case k\+ki^> k2 + ki,. In this case, 
for f < 0, two Y-shape resonances consist of solitons [2,3], [1,3], [1,2] and solitons 
[1,4],[2,4],[1,2] and for t > 0, soUtons [1,4],[1,3],[3,4] and solitons [2,3],[2,4],[3,4]. 
For t < 0, [1,4] and [2,3] solitons make an asymmetric interaction locally. This in- 
teraction generates only the height less than the height ?7i,4 of the highest asymptotic 
soliton. At f = 0, the exchange of combination of two Y-shape resonances takes 
place and four solitons interact near the origin of the xy-plane.. We can again prove 
m(0,0,0) < Ui^4 as in the resonance 2-soUton interactions. Accordingly, ?7i 4 is the 
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Figure 5: Time evolution of the inelastic 2-soliton solution obtained from Cm in ( I3.12l i with 
{ki,...,k4) = (-0.5,-0.1, 1, 1.7). Left: t = -5; center: t = 0; right: t = 5. 




Figure 6: Contour plots of the time evolution of the inelastic 2-soliton solution shown in 
Fig.|5] Left: f = -13; center: t = 0; right: f = 10. 

maximum height in this interaction. For + ^4 < ^2 + ^3 > the combinations of Y- 
shape resonances change and an asymmetric interaction of [1,4] and [2,3] solitons 
takes place for t > 0. However, the result for the maximum height is the same before. 
In the case ki + = k2 + k^, [1,4] and [2,3] solitons are parallel and make an over- 
taking interaction. At ? = 0, these two soliton may coalesce into one peak, the height 
of which is ?7i,4 — ?72,3 as described in the subsection on asymmetric 2-soliton inter- 
actions, and at the same time the exchange of combination of two Y-shape resonances 
takes place. So, in this case the maximum height is also U\^^. 

Type IV. In this case Ci,4 = 0. The incoming solitons are [1,3] and [2,4], the out- 
going solitons are [1,4] and [2,3]. Such a solution is a time-reversal version of the 
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inelastic solution of type 111: that is, uiY{x,y,t) = um^—x, —y, —t). 



3.3. Generation of large-amplitude waves 

None of the three types of elastic 2-soliton solutions describes the generation of large- 
amplitude waves from the interaction of lower-amplitude ones, since the interaction 
pattern of ordinary and asymmetric solutions is stationary, and the intermediate soli- 
ton [1,4] in resonant solutions is present at all times. On the other hand, inelastic so- 
lutions of type 1 do have the effect of an amplification of the maximum wave height. 
Moreover, when this kind of solution is embedded in a larger soliton complex, fur- 
ther increases of the wave height may result from the interaction of the interaction 
arm with the other solitons in the complex. As an example. Fig. |7] shows the large- 
amplitude wave produced by the interaction generated by the 3x6 coefficient matrix 

'^l 1 1 0^ 



C 



(3.13) 



110 
yO 1 ly 
As evident from Fig. |7J the interaction among the solitons results in the temporary 
generation of an extreme wave whose height exceeds four times that of the highest 
asymptotic soliton. 

This discussion suggests the following physical mechanism of generation of ex- 
treme waves: (i) several solitons are generated by external sources; (ii) two of those 
solitons generate a large-height interaction arm, as in the case of ordinary 2-soliton 
solutions, or that of inelastic solutions of type 1; (iii) this interaction arm interacts 
with one of the other solitons or interaction arms, in which case the wave height can 
be many times higher that of each asymptotic soliton; (iv) after the interaction, the 
wave amplitude decreases. 
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Figure 7: A solution of KP generating an extreme wave, 
as produced by the coefficient matrix in ( 13.131 ) with 
{ku-.-M) = (-2,0,2,4,4.01,8), 0i;o = ...05;O =0 and 
06;O = 200. Above, left: t = -0.4; above, center: t = 0.7; 
above, right: t — I. The plot to the right shows the maxi- 
mum wave height as a function of time. 
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3.4. A method to predict maximum interaction amplitudes from asymptotic data 

The above results yield a method to predict the possible maximum amplitude of a 
soliton interaction based only on information about the asymptotic line solitons. Con- 
sider first the case in which there are two line solitons as j — > ±00. Let {a^ ,d^) and 
{aj ,d2) be the amplitudes and directions of the incoming hne solitons, and {af ,d^) 
and (aj,<i^) be those of the outgoing line solitons, both sorted in order of increas- 
ing values of d. If d^ = d^ or d2 = d2 , the incoming and outgoing line solitons 
coincide, implying that one has an elastic soliton interaction. Otherwise one has an 
inelastic soliton interaction. For an elastic soliton interaction, the interaction ampli- 
tude is computed as follows. Let K„^± = {d^ ±<3„)/2 for « = 1,2. [The superscript 
± in the soUton parameters a„ and J„ can of course be dropped for elastic solutions.] 
Then define the phase parameters . . . ,^4 to be the values k"i -t, K"2,± rearranged in 
order of increasing size, such that /ci < ^2 < ^3 < ^4- In this way each soliton is 
uniquely identified by the index pair [injn] that labels the position of k:„ _ and k:„ + 
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(respectively) in the list /ci , . . . ,^4. The type of index overlap determines the type of 
soliton interaction |l2l|7l[T8l, and the interaction amplitude is then obtained from the 
calculations described earlier. A similar method applies in the case of an inelastic 
interaction. In this case, however, one needs the asymptotic data both as j — > —00 and 
as J —> 00 in order to uniquely identify the type of interaction. 

If the number of solitons is greater than two, one selects any two neighboring 
solitons and perform the above procedure to compute the possible maximum ampli- 
tude. One then repeats these steps for all possible pairwise combinations of solitons. 
Note, however, that, unlike the case of two solitons, this procedure does not yield a 
precise estimate, for two reasons: (i) whether or not the theoretical maximum am- 
plitude in each pairwise interaction is reahzed depends on the details of the soliton 
configuration; (ii) these high-amplitude intermediate solitons resulting from pairwise 
interactions can in some cases interact among themselves producing solitons of even 
higher amplitude. Again, whether or not this happens depends on the details of the 
soliton configuration. A true upper bound can be obtained: ?7max = ^(^max — ^min)^- 
This theoretical maximum, however, is realized only in a small number of soliton 
interactions, as should already be evident from the case of 2-soliton solutions. 

4. Numerical simulations 

We now describe numerical simulations of multi-soliton interactions of the KPII 
equation. The numerical simulation of multi-soliton solutions is particularly impor- 
tant, since at present no analytical methods exist to investigate the stability of such 
solutions using either the inverse scattering transform or other techniques. 
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4.1. A computational method for line-soliton solutions 

The numerical integration of the KP equation poses a number of challenges (e.g., see 
IITtI and references therein). In particular, when simulating soliton solutions one must 
take into account that line solitons are not localized objects, but they extend through 
the boundaries of any finite computational window. The approach we used here is 
based on the one in ll35l[36l[37l[38ll39l . but with different boundary conditions. For 
the x-direction, we set our computational window to be wide enough that any initial 
solitary waves are far away from boundary. This allows us to use periodic boundary 
conditions and to compute x-derivatives with spectral methods. For the j-direction 
we employ the windowing method |30|, which has its roots in signal processing, and 
where the windowing operation allows the spectral analysis of non-periodic signals. 
We use the following window function: W{y) = 10^'^"l^^'/^^^l " , where L is the length 
of the computational window in the j-direction, and a and n are parameters. Here we 
set a = 1 . 1 1 1 and n = 21. We then transform the solution as follows: 

q{x,y,t) =W{y)u{x,y,t), Q{x,y,t) =W{y)u^{x,y,t) . (4.1) 

Substituting this into the KPII equation, we obtain 

{-4qt+ 3Q^ + q,:,,)^ + 3qyy - 6WyUy + 3WyyU = . (4.2) 

All terms in this equation vanish at the boundaries in the 3'-direction. This makes 
it possible to apply pseudospectral methods to compute j-derivatives. We then in- 
tegrate (14.21 ) in time in the Fourier domain using Crank-Nicholson differencing and 
an iterative method. Once q{x,y,t) is obtained, u{x,y,t) is recovered from (14.11 ). But 
note that the formula for u{x,y,t) becomes ill-conditioned near the boundaries in the 
j-direction, where W{y) tends to zero. Near these boundaries, we thus correct the 
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solution using information about the soliton behavior. All the simulations were per- 
formed on a grid with 8192 x 1024 points, Ax = Ay = OA and At = 0.005. Figures 
ISUTTlbelow show the resulting field q{x,y,t). Note that, to make the interactions more 
evident, only a small portion of the computational domain is often shown. 

4.2. Numerical simulations ofline-soliton interactions 

Figure [8] shows the field q{x,y,t) corresponding to the numerical time evolution of 
an initial condition (IC) consisting of an exact fully resonant 3-soliton solution. As 
evident from the figure, the numerical solution accurately reproduces the web struc- 
ture observed in the exact solution. This result confirms that the numerical method 
described above can indeed effectively simulate multi-soliton interactions, and at the 
same time provides a first indication that such solutions are stable. 

Next we describe the time evolution of ICs consisting of a linear superposition 
of two line soli tons. Figures |9] and [TO] show the cases where the amplitudes and 
directions were chosen so as to correspond respectively to an ordinary and resonant 
interaction. We emphasize that, in both cases, the initial state is not an exact two- 
soliton solution. Indeed, in both cases the numerical solution shows the presence 
of radiative component of small amplitude. Nonetheless, the results do provide a 
further check of the stability of 2-soliton interactions. Note in particular that the 
characteristic "box" of the resonant solution is generated numerically. Similar results 
were obtained by numerically computing the time evolution of an IC consisting of a 
linear superposition of two line solitons corresponding to an asymmetric interaction 
and of three line solitons corresponding to a fully resonant solution. Finally, Fig. [TT] 
shows the time evolution of an IC corresponding to an inelastic interaction. 

Importantly, when the IC is not an exact 2-soliton solution, the numerical solu- 
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tions show that the height of the interaction arm tends monotonically in time to the 
value for the corresponding exact soliton as obtained in section |3]] For example, in 
Fig. |9] the height of interaction arm increases monotonically in time from its initial 
value of 2 [owing to an IC that is just a linear superposition of two line soli tons], ap- 
proaching asymptotically the value corresponding to ordinary 2-soliton interactions. 
Conversely, for an asymmetric solution the height of the interaction arm decreases in 
time, again approaching asymptotically the value for asymmetric 2-soliton interac- 
tions. These results extend the validity and usefulness of the analysis in section [3l 

The above results suggest that multi-soliton solutions of KPII are robust and sta- 
ble, and, moreover, that even when the solution contains non-solitonic components, 
only the information about the asymptotics line solitons contributes to determine the 
intermediate interactions of the line solitons. We also suspect that, as the radiative 
components disperse, asymptotically in time the solution will be closely approxi- 
mated by an exact soliton solution, similarly to the case of (l-i-l)-dimensional soliton 
equations. Of course all of these conjectures must be carefully tested and validated 
with extensive numerical simulations, which are beyond the scope of this work. 

5. Conclusions 

We have studied the amplitude of soliton interactions of the KPII equation, and we 
have discussed a possible mechanism for the generation of large-amplitude waves. 
Ordinary A/^-soliton interactions with N >3 can also briefly produce large amplitude 
waves if all the solitons intersect simultaneously at the same point in the xy-plane. 
Note however that the event of all solitons intersecting simultaneously at a single 
point can be considered to be statistically speaking unlikely in a multi-soliton com- 
plex, with a likelihood decreasing as the number of solitons increases. In this sense, 
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therefore, the mechanism described in this paper, involving inelastic 2-soliton solu- 
tions (possibly embedded in a larger soliton complex, as in Fig. |7]l represents the most 
likely way to generate large-amplitude waves. 

We also proposed a method to determine the maximum amplitude resulting from 
the interaction of two line solitons. The calculation of the maximum amplitude is 
based on the framework of exact line-soliton solutions, but it may also be useful for 
solutions where a non-solitonic component is present, if line-soliton solutions of KPII 
are indeed proven to be stable, since then, even if one starts from an initial state that is 
not an exact soliton solution, the radiative portions of the solutions will disperse away, 
and asymptotically in time one will approach a state consisting of an exact soliton 
solution. An interesting open problem will be to develop an algorithm to compute the 
actual maximum amplitude generated by any multi-soliton configuration. 

Finally, we implemented an algorithm to numerically integrate solutions of the 
KPII equation containing multi-soliton complexes, and we discussed the results of 
numerical simulations. These results show the robustness of all types of line-soliton 
solutions of KPII, including those exhibiting web-like structure. We also confirmed 
numerically that multi-soliton interactions can generate large amplitude waves, and 
that an initial state that is not an exact solution eventually converges to an exact multi- 
soliton solution. We note that resonant 2-soliton solutions with a hole have also been 
found in other discrete and continuous (2-1-1 )-dimensional soliton equations ifTSlfTSl 
[T9l|2Tl[35l, both in exact solutions and in numerical simulations, indicating that they 
are a fundamental and robust structure of (2-i-l)-dimensional soliton equations. 
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Figure 8: Numerical time evolution of the fuUy resonant 3-soliton interaction; {ki,. . . jk^) = 
(-2.5,-1.5,-0.5,0.5, 1.5,2.5). Left: t = 0; right: t = 10. 




Figure 9: Numerical time evolution of a linear superposition of Une soUtons generating an 
ordinary 2-sohton interaction; {ki,... ,^4) = (-1,-0.001,0, 1). Left: t = 0; right: t = 50. 
The IC is not an exact 2-soUton solution, and some dispersive waves are generated, but the 
results suggest the stability of the sohton interactions. 




Figure 10: Numerical time evolution of a Ij^ar superposition of hne solitons generating 
a resonant 2-soliton interaction; {ki,... ,^4) = (-2,-1, 1,2). Left: ^ = 0; right: t = 12.5. 
Again, dispersive waves are generated, but the interaction appears to be stable. 



Figure 11: Numerical time evolution of an inelastic 2-soliton interaction. Left: t = 0; center: 
t = 2.5; right: t = 5. 
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